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Abstract 

We extend a recent construction for an integrable model describing 
Josephson tunneling between identical BCS systems to the case where 
the BCS systems have different single particle energy levels. The exact 
solution of this generalized model is obtained through the Bethe ansatz. 



1 Introduction. 

The experimental work of Ralph, Black and Tinkham|l], ^ on the discrete energy 
spectrum in small metallic aluminium grains has generated substantial interest 
in understanding the nature of superconducting correlations at the nano-scale 
level. Their results indicate significant parity effects due to the number of 
electrons in the system. For grains with an odd number of electrons, the gap in 
the energy spectrum reduces with the size of the system, in contrast to the case 
of a grain with an even number of electrons, where a gap larger than the single 
electron energy levels persists. In the latter case the gap can be closed by a 
strong applied magnetic field. The conclusion drawn from these results is that 
pairing interactions are prominent in these nano-scale systems. For a grain with 
an odd number of electrons there will always be at least one unpaired electron, 
so it is not necessary to break a Cooper pair in order to create an excited state. 
For a grain with an even number of electrons, all excited states have a least one 
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broken Cooper pair, resulting in a gap in the spectrum. In the prescence of a 
strongly applied magnetic field, it is energetically more favourable for a grain 
with an even number of electrons to have broken pairs, and hence in this case 
there are excitations which show no gap in the spectrum. 

The physical properties of a small metallic grain are described by the reduced 
BCS Hamiltonian[§ § 

c c 

Hbcs = e i n i ~ 9 4+4- c J- c i+- W 

Above, j — 1,...,C labels a shell of doubly degenerate single particle energy 
levels with energies ej and rij is the fermion number operator for level j. The 
operators Cj±, cj ± are the annihilation and creation operators for the fermions 
at level j. The labels ± refer to time reversed states. 

One of the prominent features of the Hamiltonian ([!]) is the blocking effect. 
For any unpaired electron at level j the action of the pairing interaction is zero 
since only paired electrons are scattered. This means that the Hilbert space 
can be decoupled into a product of paired and unpaired electron states in which 
the action of the Hamiltonian on the subspace for the unpaired electrons is 
automatically diagonal in the natural basis. In view of the blocking effect, it 
is convenient to introduce hard-core boson operators bj — Cj-Cj+, b\ = cj + cj_ 
which satisfy the relations 

(&]) 2 = 0, [b 3 ,b{]=5 3k (l-2b%) 

[b j ,h] = [b],bl} = 

on the subspace excluding single particle states. In this setting the hard-core 
boson operators realise the su(2) algebra in the pseudo-spin reprepresentation, 
which will be utilized below. 

The original approach of Bardeen, Cooper and Schrieffer|| to describe the 
phenomenon of superconductivity was to employ a mean field theory using a 
variational wavefunction for the ground state which has an undetermined num- 
ber of electrons. The expectation value for the number operator is then fixed by 
means of a chemical potential term /x. One of the predictions of the BCS theory 
is that the number of Cooper pairs in the ground state of the system is given by 
the ratio A/d where A is the BCS "bulk gap" and d is the mean level spacing 
for the single electron eigenstates. For nano-scale systems, this ratio is of the 
order of unity, in seeming contradiction with the experimental results discussed 
above. The explanation for this is that the mean-field approach is inappropriate 
for nano-scale systems due to large superconducting fluctuations. 

As an alternative to the BCS mean field approach, one can appeal to the 
exact solution of the Hamiltonian ([[]) derived by Richardson and ShermanQ. It 
has also been shown by Cambiaggio, Rivas and SaracenoQ that ([j]) is integrable 
in the sense that there exists a set of mutually commutative operators which 
commute with the Hamiltonian. These features have recently been shown to 
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be a consequence of the fact that the model can be derived in the context of 
the Quantum Inverse Scattering Method (QISM) using a solution of the Yang- 
Baxter equation associated with the Lie algebra su(2)|| 9. One of the aims 
of the present work is to extend this approach for application to generalised 
models. As a specific example, we will show that a model for strong Josephson 
coupling between two BCS systems falls into this class. 

Recall first that electron pairing interactions manifest themselves in macro- 
scopic systems via three well known phenomena: 

• supercurrents 

• Meissner effect 

• Josephson effect 

As noted by von Delft the notion of a supercurrent in a nano-scale system 
is inapplicable because the mean fee path of an electron is comparable to the 
system size. Likewise, the penetration depth of an applied magnetic field is 
comparable to the system size, which prohibits any Meissner effect. 

Josephson [[fl| put forth a proposal for the tunneling of electron pairs between 
superconductors separated by an insulating barrier. A theory was derived to 
describe weak coupling between two superconductors treated at the mean field 
level in the grand-canonical ensemble. A remarkable prediction of the theory 
was that it is possible for a direct current to flow across the insulator for the 
case of zero applied voltage, whereas a constant voltage across the insulator 
produces an alternating current. The essential features of the theory stem from 
the phase difference between the superconductors, which is well defined since the 
variational wavefunctions for the superconductors have undetermined particle 
numbers. 

For the case of nano-scale systems, the above predictions are again invalid 
due to the finite particle numbers for each system, giving rise to phase uncer- 
tainty. However, if we are to consider strong coupling where individual particle 
numbers are not conserved, only total particle number, it is appropriate to study 
the effective Hamiltonian 

c 

H = H B cs(l)+H BCS (2)~sjY,(b](l)b k {2) + b]{2)b k {l)y (2) 

where ej is the Josephson coupling energy, for the purpose of investigating the 
nature of pair tunneling at the nano-scale level. In a previous workjl^] it was 
shown that the above Hamiltonian is integrable for ej = g for the case when 
Hbcs{^): Hbcs(2) have identical single electron energy levels. Below we will 
extend this construction to the case where Hbcs{^)> Hbcs{2) describe non- 
identical systems. 
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2 A universal integrable system. 

First we introduce the Lie algebra su(2) with generators S + , S~, S z satisfying 
the commutation relations 

[S z , 5 ± ] = ±S ± , [S+, S~] = 2S Z . (3) 

The Casimir invariant, which commutes with each element of the algebra, has 
the form 

C = S + S~ + S~S + + 2{S z f . 

Associated with the su(2) algebra there is a solution of the Yang-Baxter equation 
in End+®EndV<g)su(2), where V denotes a two-dimensional vector space. This 
solution reads ]T^| 

Ru{u - v)Li(u)L 2 (v) = L 2 (v)Li(u)Ri 2 (u - v) 

with 

2 

R(u) = JSJ+^e"®^, 

m,n 

L{u) = 7(g)7+ ^{e\®S z -el®S z +e\®S-+el®S + ) 

where {e™} are 2x2 matrices with 1 in the (m, n) entry and zeroes eleswhere. 
Above, 7 is the identity operator and rj is a scaling parameter for the rapidity 
variable u which plays an important role in the subsequent analysis. With this 
solution we construct the transfer matrix 

t(u) = tr (G L 0C (u - e c )--L m {u - e x )) (4) 

which is an element of the £-fold tensor algebra of su(2). Above, tro denotes the 
trace taken over the auxiliary space and G = exp(ary(r) with a = diag(l, — 1). 
A consequence of the Yang-Baxter equation is that [t(u), t(v)] = for all values 
of the parameters u and v, and independent of the representations of su(2) in 
the tensor algebra. Defining 

Tj = lim o^tiu) 

for j = 1, 2, £, we may write in the quasi-classical limit Tj = tj + o(j]) and it 
follows that [Tj, 77-] = 0, Vj, k. Explicitly, these operators read 

r 3 =2aS z + J2-^- (5) 

W 3 k 

with 9 = S+ ® S- + S- ® S+ + 2S* Z ® S 12 . 
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We define a Hamiltonian through 



3 j-k j j 



E 2 ^j-^XXs+ (7) 



J-k 



The Hamiltonian is universally integrable since it is clear that [H, Tj] — 0, Vj 
irrespective of the realizations of the su(2) algebra in the tensor algebra. 
Realizing the su(2) generators through the hard-core bosons; viz 

St = bj, Sj = b\, S*= l -{I-n 3 ) (8) 

one obtains (Q) (up to a constant) with g = 1/a as shown by Zhou et al.[[| and 
von Delft and PoghossianQ. 

We now turn to applying (f7j) for the study of two coupled BCS systems. To 
accomodate this, it is convenient to first consider three index sets Po, Pi, P2 
such that individually the BCS Hamiltonians are expressible 

c c 

H B cs(i)= e i n i-9 E b k b i- 

je(PoUPi) i : fcG(-P up I ) 

If the single particle energy Cj is common to both systems, then j £ Pq. Hence 
it is meant to be understood that ej ^ e*. ^ e;V j G Pi, k G P2, I G P3. In the 
case that j G Po, the local su(2) operators are described by the tensor product 
of two pseudo-spin realisations acting on the four-dimensional tensor product 
space. We can now realise (^) in terms of the hard-core boson representation 

# = &,•(<), Sr=b](i), S* = -(I-nj(i)) 
for j G Pj, i = 1,2 whereas for j G Po we take the tensor product representation 

S+ = 6,(1) +6,(2) 
S i = b}(l)+b](2) 
S? = 7-^(^(1) +^(2)). 

Under this representation of (Q) we obtain (||) with sj = g = 1/a, establishing 
integrability at this value of the Josephson coupling energy. For the case when 
the index sets Pi, P2 are both empty, i.e., the two BCS systems are identical, 
this result was previously shown by Links et al. fllQf . 
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3 The exact solution. 



In addition to proving integrability for sj = g, we can also obtain the exact 
solution from the Bethe ansatz. Below we will derive the energy eigenvalues for 
the Hamiltonian (Q) in a very general context, which includes those of (||) with 
Ej = g as a particular case. 

For each index k in the tensor algebra in which the transfer matrix acts, 
and accordingly in (Q), suppose that we represent the su(2) algebra through the 
irreducible representation with spin s k - Thus {S£, , S%} act on a (2s k + 1)- 
dimensional space. Employing the standard method of the algebraic Bethe 
ansatz gives that the eigenvalues of the transfer matrix (|j) take the form 

A(«) = eM*v)f\ U ~ ek + 71Sk f\ U ~ W3 ~ T1 

J - l - L U — Ck ±A - U — Wn 

k 3 



C M 



exp 



U — €k — r]s k yj u - Wj + r\ 



U — E k U — Wj 

k J ■> 



Above, the parameters Wj are required to satisfy the Bethe ansatz equations 



M 



Wl- Cfe + T]S k __ j-rWl-Wj+Tj 



exp(2„,)n ::"; r =-n 



, wi - e k - rjs k -V- 1 - Wi - Wj -T] 

K J 

The eigenvalues of the conserved operators (||) are obtained through the 
appropriate terms in the expansion of the transfer matrix eigenvalues in the 
parameter rj. This yields the following result for the eigenvalues Xj of Tj 

such that the parameters Vj satisfy the coupled algebraic equations 

2a + Y^-=jr^—. (10) 

Through (|) we can now determine the energy eigenvalues of (Q) . It is useful to 
note the following identities 

M M C 

2 «E^ 2 EE^ = m{m — i) 

3 j k 3 

M C 

aM+yy-^- = o 

3 k 3 

M C M C C 

EE^-EE^ - M E si - 

3 k J 3 k J k 
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Employing the above it is deduced that 



3 j 

C C C C CM 

X>Aj - 2a^e jSj +^^ SjSfc -2M^ Sfc -2a^^+Af(Af-l) 

3 3 3 k^j k j 

which, combined with the eigenvalues 2sj(sj + 1) for the Casimir invariants Cj, 
yields the energy eigenvalues 

M C 

£ = 2E^- 2 E s ^- (ii) 



From the above expression we see that the quasi-particle excitation energies are 
given by twice the Bethe ansatz roots {vj} of (|l0|). 

In order to specialise this result to at integrable coupling, it is useful to 
first make the following observation. For j £ Pq, in which case the su(2) algebra 
is realised via the tensor product of two hard-core boson representations, it is 
well known that the representation space is completely reducible into triplet 
states and a singlet state. Note however, that for the singlet state the su(2) 
generators act trivially, and hence this state is blocked from scattering in analogy 
with the blocking of single particle states discussed in the introduction. Hence 
the su(2) algebra will only act non-trivially on the triplet states. In specialising 
( [TO 11) to the case of (^), we need only to set Sj = 1/2 for j e Pi U P2 and 



= 1 for j G P . 



4 Conclusion 

We have displayed the existence of a general class of integrable systems which 
includes the reduced BCS Hamiltonian and a model for strong Josephson tun- 
neling between two reduced BCS systems. By deriving the models through 
the QISM we have also determined the exact solution via the Bethe ansatz. 
A further application of this approach is the computation of form factors and 
correlation functions^, [T?| . 
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